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Abstract 

We examine the role of matter in an oscillating background with a warped, 
compact extra dimension. This background is compatible with an S 1 /Z2 orb- 
ifold structure which allows chiral fermions to be included in the scenario. 
When the background oscillates rapidly, the leading coupling of these oscilla- 
tions is to gauge fields rather than fermions. If the decay of these oscillations 
were to occur today, it could provide an alternative mechanism for generating 
the ultra high energy cosmic rays. 
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1. Introduction. 

A remarkable property of actions that generalize the Einstein-Hilbert action for gravity 
is that they admit metrics with a periodic dependence on one or more of the coordinates. 
The simplest such example occurs for an action composed of a general set of curvature 
invariants, with up to four derivatives of the metric, and scalar field. In 4 + 1 dimensions 
a metric periodic in one of the spatial dimensions provides a naturally compact space — 
by identifying the compactification length with the period — without any singularities or 
discontinuities in the metric [0. Moreover, these periodic solutions exist without the need 
for finely tuning any of the parameters in the action. The coefficients of the four-derivative 
terms determine an essentially unique compactification size. 

Similarly, we can search for metrics that oscillate both in the fifth dimension and 
in time [0. Although these metrics do not describe the most general isotropic vacuum 
solutions for the action, they could have intriguing consequences for the early and the late 
evolution of the universe. In the very early universe, we shall see that the rapid decay of 
the oscillating background into gauge fields provides an example of how a universe with a 
Planck-scale time dependence can relax into one that can admit a realistic, slow evolution. 
Moreover, if some slowly relaxing, but rapid, oscillation persisted or arose recently, then 
it could contribute a flux of ultra high energy cosmic rays. To be successful, the theory 
would still require a mechanism to maintain these oscillations at a satisfactory rate of 
decay today. 

In this article we shall explore the effect of including fermion and gauge fields in 
these oscillating backgrounds. Even for a purely static background, introducing chiral 
fermions requires that the compact extra dimension should have the topology of an S 1 /^ 
orbifold. The metric and the scalar field that supports the compact geometry must be 
respectively even and odd under this Z2. When the metric additionally oscillates rapidly 
in time, the coupling of these fields with gravity leads to a potential decay channel for 
the oscillations. We shall show by expanding a fermion and an Abelian gauge field in 
Kaluza-Klein modes that to leading order the fermion zero mass mode does not couple to 
the oscillating component of the metric. The preferred channel for the relaxation of the 
metric is therefore into gauge fields. 

In the following section we review the features of oscillating metrics in 4+1 dimensions. 
Section 3 describes how to place the theory on an orbifold to allow for chiral fermions in 
a warped, static background. Section 4 examines the coupling of the time-dependent 
components of the metric with the Kaluza-Klein modes of a fermion and an Abelian gauge 
field placed in an oscillating, compact background. Finally, in section 5, we study the 
phenomenological signature in ultra high energy cosmic rays of a very small, but rapid, 
oscillation today. 



2. Oscillating metrics. 

At distances approaching the Planck length, the usual Einstein-Hilbert action should 
be supplemented by higher derivative curvature invariants. For example, an action with 
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up to four derivatives of the metric and a scalar field has the formEl 

S = Ml J d 4 xdy [-2 A + R + a R 2 + b R ab R ab + c R abcd R abcd ] 
+ Ml j d 4 xdy [- 1 V a 0V°0 + A£ J +■■■. 



(2.1) 



M5 and A denote the Planck mass and cosmological constant respectively. This metric 
admits smooth, non-singular metrics that are compact in the fifth dimension, 

ds 2 = g ab dx a dx b = e A{y) dx»dx u + dy 2 (2.2) 

where the exponent A(y) is a periodic function, when is either an interaction, such 
as (V a (pV a (j)) 2 HJ, or a Casimir term with differing values its components in the large 
ordinary x^ and the compact y directions. In either case, the existence of such solutions 
does not require finely tuning any of the parameters in the action. The effects of the R 2 
terms can be conveniently parameterized by 

^ = 160 + 56 + 40, A = 5a + b + \c and v = 3a + b + c. (2.3) 

A and u, in particular, represent the coefficients of a Gauss-Bonnet and a squared Weyl 
tensor respectively. 

While the static metrics are adequate for studying a theory with a flat 3+1 dimensional 
long distance limit, models that are to include a realistic cosmology should evolve in time 
as well, 

rfs 2 = _ e A(y) dt 2 + e A(y) e B(t) g.. ^i^j + g C(t) dy 2. ^ 

here we have still assumed an isotropy in the three large spatial dimensions but have 
allowed a different evolution in the compact dimension. 

We shall principally consider a universe that oscillates rapidly so that the time de- 
pendence and the y-dependence, which fixes the size of the compact dimension, are on 
a similar footing. The fourth-order, two-variable differential equations that result from 
varying the action (2.1) are more difficult to solve than the static case (2.2), which can 
be solved numerically, but fortunately it is possible to build a solution order by order in a 
small amplitude expansion. To first order in e y , et ^ 1, we find 

Mv) = £y cos(uj y y) + ■■■ 

(2.5) 

B{t) = et cos(u>tt) + ■ ■ ■ C(t) = Set cos(uJtt) + ■ ■ ■ 

with 



U y = \lj- and Ut= V3/-16, (2 - 6) 



1 Our convention for the signature of the metric is ( — ,+ ,+ ,+ ,+) while the Riemann curvature 
tensor is defined by -R a bcd = d d T% c -d c r% d +Tl d Tt c -r% c rt d . Here the indices a, b, . . . = 0, 1, 2, 3, y 
range over all coordinates while fj,, v ... = 0, 1, 2, 3 label the ordinary space-time coordinates. 
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so that a periodic solution exists when fi < and v < j^fx — no fine-tuning is needed. In 
(2.6) we see explicitly that the higher derivative terms in the action determine the size 
of the extra dimension and the oscillation frequency. Proceeding to the next order in the 
small e y , et expansion introduces some corrections to each of A(y), B(t) and C(t), which 
depend on both t and y, and more importantly relates the small amplitudes to the size of 
the cosmological constant: 

A=fM4+^ e *)- ( 2 - 7 ) 

The relative size of the e y and et in this solution is undetermined. This relation allows us to 
state more precisely the regime in which the small amplitude expansion exists. When p and 
v have a natural size, 0{M§ 2 ), then a small amplitude translates to a small cosmological 
constant, M5A -1 / 2 ^> 1. Note however that the existence of periodic solutions — beyond 
the small amplitude regime — does not require a small cosmological constant as was shown 
numerically in [[TJ. 

In the large 3 + 1 dimensions, once we have substituted the leading effects from the 
y-dependent scalar field, we discover that the density p and the pressure p for the non- 
compact spatial dimensions are 

p = \u\e\ + \4> 2 + ■ ■ ■ « f u 2 e 2 > 

(2.8) 

p = -^u 2 e 2 + U 2 + ---^0 



so that \ui 2 e 2 resembles an effective cosmological constant. The fundamental cosmological 
constant that appears in the action effects the warping of the extra dimension, as the 
pressure in that dimension indicates, p y = — 2 A + ^u 2 e 2 . 

We find that the small amplitude expansion described above cannot be extended to 
the case of a de Sitter expansion B(t) = kt + ■ ■ ■ that coexists with the oscillations, if the 
rate of expansion is to be comparable to the leading oscillating terms in the metric, i.e. 
k ~ 0(M$et). In general, the dynamics that links the short scale oscillations with the 
large scale evolution of the universe is complicated and deserves further study. 

If the amplitude of the oscillations starts with a value large compared to any expansion 
rate, the amplitude must eventually decrease sufficiently so as to allow a realistic cosmology. 
Given that the vacuum energy density in (2.8) is positive, a natural mechanism for its 
relaxation is through the decay of the rapid oscillations into energetic particles, analogous 
to the decay of the infiaton in inflation. In section 4, we shall study how fermions and 
gauge fields couple to the oscillating functions, B(t) and C(£), thus providing an explicit 
decay mechanism. 

Note that the presence of the extra dimension is necessary for this decay. In 3 + 1 
dimensions the metric would not contain a C(t) term and would therefore be conformally 
flat; then the oscillation could not decay into conformally coupled fields. 



3. Chiral fermions. 

In order to recover the standard model fields, the theory must contain chiral fermions. 
A difficulty arises for the usual Kaluza-Klein compactification with only one extra dimen- 
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sion since the Lorentz symmetry group SO(4, 1) has only one, non-chiral, spin-| represen- 
tation. Chiral spinor representations do arise when we break the full SO(4, 1) symmetry 
group, for example by placing the theory on an S 1 /^ orbifold in the fifth dimension H]. 

We can similarly introduce an orbifold into the warped Kaluza-Klein picture with a 
static metric, (2.2). The field equations for the action (2.1) that determine A(y) and <p(y) 
do not depend explicitly on y so we are free to translate an extremum of A(y) to y = 0. If 
the period of A(y) is y c then we can compactify the space by restricting y e [— hy c , hy c ] and 
identifying the endpoints. The solutions for A(y) found numerically in jl| and analytically 
in a small amplitude expansion in || are manifestly even under y — > — y so the background 
metric respects an SO (3, 1) x Z2 symmetry. 

This Z2 invariance allows us to introduce an orbifold geometry in the extra dimension 
by identifying y with —y. In order to define a theory consistently on this orbifold, the fields 
must be odd or even under this discrete Z2. As we have seen, the background metric is 
even and this orbifold geometry will further constrain the allowed gravitational excitations 
of this background. 

The scalar field, in contrast, must be odd, <fi(—y) = —<p(y). The reason is that the 
field equations for (2.1) relate <fi'(y) to an expression that depends on A(y) only through 
terms with even numbers of derivatives. Therefore, (f)'(y) oscillates with the same period 
as A(y) and is also even under <p'(y) = <p'(—y). For the solutions found in 0, <p'(y) 
is everywhere positive so that after integrating we obtain a (j){y) that increases monotoni- 
cally. Choosing the constant of integration so that (f)(0) = produces a (f)(y) that is odd. 
To accommodate the boundary values, we must further impose that <p(y) itself assumes 
values only in a compact space by identifying cf>(y c /2) and (f)(—y c /2). In figure 1, we show 
explicitly a example with this geometry. The form for 4>(y) was found numerically for 
AC^ = k(V a (pV a 4>) 2 but with the parameters chosen arbitrarily within the region of the 
{A, \i, A, /c}-space with periodic solutions in y. 




Figure 1. The behavior of <f>(y) on the orbifold for A = 1, fj, = 0.1, A = 0, and AC^ = 
a (f>y a <t>f ■ The initial condition is A"(0) = 23.77364592. Here y 6 [-0.6578,0.6578] and 
4> 6 [—0.744,0.744]. The endpoints of each dimension are identified so that both are compact. 

In order to be compatible with this Z2 orbifold structure, a separate scalar field should 
be included if we wish also to allow the rapid time oscillations. This requirement follows 
from the leading order behavior, <fi 2 = 3u) 2 6 2 @; upon integration we obtain a contribution 
that is even under the Z2 symmetry which is incompatible with the symmetry of the scalar 
field that supports the compact extra dimension. 



4 



We now obtain chiral fermions through the standard construction 0] [BJ • If we choose 
the following boundary conditions on a five dimensional fermion, 

^ A ,- !/ )=7^(x A ,»/), (3.1) 

then expanding in a tower of Kaluza-Klein modes, we have a chiral zero mass mode, 

^\x\y)=0 , ^\x\y) = e- A ^^\x x ), (3.2) 

as well as a series of paired massive left and right modes ip^ R (x x ,y) for n > 0. Unlike 
the flat Kaluza-Klein expansion, the zero mass mode does depend on the fifth coordinate 
through the e~ A ^ factor. 



4. Fermions and gauge fields in an oscillating background. 

The presence of ordinary fermion or gauge fields affects the evolution of universe with 
rapid time oscillations. The interaction of such fields with the oscillating metric offers a 
route for the relaxation these oscillations. However, we shall see that the zero mass mode of 
the Kaluza-Klein expansion of a massless five dimensional fermion does not couple directly 
to the oscillating components of the metric, to leading order. A decay into gauge fields (or 
scalar fields) becomes the dominant channel for the relaxation of the metric. 

Returning to a metric with a dependence both on time and the extra dimension (2.4), 
when B(t) ^ C(t) this metric is not conformally flat. This feature becomes more apparent 
if we define new coordinates 

V (t) = J e'^'Ut' and r(y) = ^ e ~W> dy', (4.1) 

in terms of which the metric (2.4) becomes 

ds 2 = e A{r) e B ^ \-drf + S tJ dx l dx J + e c ^- B ^ dr 2 } , (4.2) 

where A(r), B(rj) and C(rj) are understood to be A(y(r)), B(t(rj)) and C(t(rj)). 

We now introduce a fermion field tp and an Abelian gauge field A a through the actioni 

S^a = J d"xdr y/=jj {e^T A (D a - igA a )^ - \g ab g cd F ac F bd ) , (4.3) 

where F a i is the gauge field strength. We work in the limit in which these fields do not 
significantly affect the background geometry. In this action, we have not assumed an 



2 Because of our choice for the metric's signature, we have not included the usual factor of i 
in the fermion action. The 7 matrices are defined by T A = {^7^,75}, where 7^ and 75 are the 
standard 7 matrices, so that {T A ,T B } = 2r) AB for r] AB = diag[-l, 1, 1, 1, 1]. Refer to § or § for 
further details on fermions in a curved background. 
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orbifold structure in the extra dimension, but as in the previous section one can be readily 
introduced. The e£ is a vierbein which connects the curved coordinates with a locally 
Lorentzian frame: 

e£e£g ab = n AB . (4.4) 
The covariant derivative of a fermion, 

D a = d a + \u aBC cj BC , (4.5) 

involves the spin-connection defined by 



uJaAB = \e\ [d a eBb - d b e Ba ] ~ \e B [d a eAb ~ d b e Aa ] 
~ \ e A e B [d c eDd - d d e Dc ] e D a 



(4.6) 



and a BC = \{T B ,T C ]. 

The contributions from the spin connection are canceled by introducing a rescaled 
fermion field defined by 

jj(x x , r) = e- A ^e--* BM e--* c ^q(x\ r). (4.7) 
The fermion-gauge field action then becomes 

S AA = [ d 4 xdr |z^7 M (<9 M - igA^ + e^ B - c ^ 7 5 (d r - igA r )^>) 

J , (4-8) 
+ J d*xdre^ A+c ^ {-^ u V Xp F, X F up - \e^ B ~ c \^ F. r F vr } ; 

note that the first term no longer contains A(r), B{rj) or C{rf). 

4-1. The Kaluza-Klein expansion. 

To investigate the effect of the oscillating metric in the long distance limit, we expand 
the fermion and gauge fields in a tower of Kaluza-Klein modes. To simplify, we work in 
the axial gauge A r = 0. The fermions are expanded separately in left- and right-handed 
fields, 

oo 

* lA*\ r) = \(lTl^ = Yl *£i(* A )/i?i(0 ( 4 - 9 ) 
where /^^(r) satisfy 

drf ( L n) = mn/W dj^ = -m n f^ (4.10) 
and obey the following orthogonality condition, 

dr/i m) >)/i n) (r) = f ^ dr f™*{r)f£\r) = (4.11) 

J-r c /2 J~r c /2 
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where r c is the volume of the extra dimension. Analogously, we expand the gauge field 

oo 

A M (x A ,r) = ^4^(x A )^)(r), (4.12) 



n=0 

defining the masses of the modes through 



d r 



e^ A dM n A = -M%e* A hW (4.13) 



and normalizing the states through 

rr c /2 

/ dre* A{r) h {rn) {r)h {n) {r) = 5 mn . (4.14) 

J-r c /2 

The fermion-gauge interaction will induce couplings among the various Kaluza-Klein 
modes, 

GZr = 9 drf^(r)f^ R (r)h {p) (r). (4.15) 

J-r c /2 

The effective action that appears in four dimensions as a result of (4.8) is thus 

s*,a = [ ^{E^/^-^^?)^ 

m,n p 

+ e *** m v (? mn ^ - * E GT p ^y^ 

m,n p 

+ e h( B ~ c ) m n (*i n) *? + tf^tfg ) } 

n 



(4.16) 



At low energies, well below the Planck scale, only the dynamics of the massless modes 
will be important in the effective theory. The r-dependent parts of the lowest-lying fermions 
are 

/<»> = /<»> = r -v» 

along with a factor of e~ A ^ from (4.7). The couplings among the massless modes simplify 
to 

go = Gf° = GT = °- drh^(r). (4.17) 

r c J-r c /2 

The low energy effective action is thus 

S eS = J d'x {^<Vfa. - ^o4 0) )* (0) - \e* c F$F<* s >'»' + •••}. (4.18) 



In this effective action, the gravitational oscillations do not couple directly to the 
fermions, but rather to the gauge fields. In the small amplitude limit that we have assumed, 
the leading interaction is 

interaction = ~\C (t)F$ F^ ^ . (4.19) 

This coupling offers a natural channel for the decay of the oscillating component of the 
metric. Similarly, any fundamental scalar fields in the theory would also allow the decay of 
the oscillating gravitational field. The kinetic energy term for the zero mass Kaluza-Klein 
mode of a massless five dimensional scalar field has a prefactor of e B+ ^ c . 

If the fermion in (4.18) has a small realistic mass m, either through a fermion mass 
in the five dimensional action or if a massless fermion subsequently develops a dynamical 
mass, a mass term would produce a coupling between the fermion and the gravitational 
oscillations. However, this term would be suppressed by m/M^ relative to (4.19). 

5. Ultra high energy cosmic rays. 

If some slowly decaying residual oscillations were to exist today, the gauge field prod- 
ucts of this decay would provide a possible source for the ultra high energy cosmic rays. 
Here we shall only focus on demonstrating that a realistic cosmic ray spectrum can arise 
and establishing a rough limit on the rate of the decay and shall not attempt to develop 
a detailed model that produces this rate. The dominant signal would likely result from 
a decay of the oscillating background into a pair of gluons. An initial gluon of energy 
Eq ~ Mpi ~ M5 will fragment into a high multiplicity jet of particles with a wide range of 
energies. For an initial photon pair the secondaries are only those produced as the photons 
scatter from the intergalactic radio background, and thus should not provide as strong a 
constraint as gluon-gluon production. Additionally, the presence of any gauge interactions 
beyond the standard model could provide other decay channels. 

The observed flux of ultra-high energy cosmic rays, assumed here to be protons, can 
then set a bound on rate of production, or equivalently on the rate of decay of the effective 
vacuum energy density, p v3iC . We estimate the flux of protons above some energy E > by 

J^^p^E^iN,. (5.1) 

N > is the number of protons with energy E > E > in a jet produced by the initial gluon of 
energy Eq ~ Mpi. I is the attenuation length of these protons due to scattering from the 
cosmic microwave background radiation |§. When E > « 10 11 GeV then £ is approximately 
a few tens of Mpc 0. This energy is just above the expected but not seen GZK cutoff ||, 
which is exceeded here by high energy protons originating within a distance I of us. 

We shall obtain an estimate of N > from the perturbative analysis of multiparticle 
production in jets based on the modified leading log approximation (e.g. [|TD| ).eI The 

3 In this approximation the main difference at leading order between a quark or gluon jet is 
that the multiplicity of particles in the gluon jet is enhanced by a factor of | relative to the quark 
jet for the same initial jet energy. The mean particle energy of products will correspondingly be 
slightly lower in the gluon jet. In the following we have treated a gluon jet the same as a quark 
jet. 
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resulting "limiting spectrum" is known but since the initial gluon energy is so high, r = 
ln(i?o/AQCD) 3> 1, we may simplify further with a Gaussian representation. We extract 
the following results from [|10| where £ = \u.(Eq/E). 

dN ( i(£-f) 2 \ 

cxexp (5.2) 



(5.3) 



3 



3^ V 4z 




Here z = y/16N c T/b, C = a 2 /(16N c b), 3b = 11N C - 2n f , and 3a = 11N C + 2n f /N*. We 
also note that a Gaussian spectrum resembles the results from Monte Carlo simulations 
ig]. For the total multiplicity Af = J™(dN/d£) d£ we use the full "limiting spectrum" 



result, 

AT = r(S)^j J B+1 (z), (5.4) 

where B = a/b. This gives M ~ 7 x 10 5 . Assuming that 5% of the energy from each 
initial gluon emerges as protons, these results imply that N > m 3000 for E > = 10 11 GeV 
and rif = 6. This value for iV> triples when the 0(t -1 / 2 ) corrections in (5.3) are absent. 
Another indication of the sensitivity of the result to the approximation is that no single 
value of rif is actually correct over the range of energies involved, and N > ranges from 
1400 to 6000 as n f ranges from 3 to 8. This last result also shows how N> can be strongly 
affected by new physics beyond the standard model. 

The observed integrated flux of cosmic rays with E > 10 11 GeV is approximately 



2 x 10" 20 cm -2 s" 1 sr _1 0, and so (5.1) and iV> w 3000 implies that 

p vac < 3 x 10 _53 gcm- 3 s _1 . (5.5) 

Thus the current vacuum energy density of the universe, assumed to be about two-thirds 
of the critical density, has a decay rate 

^£<l -6_ (5. 6 ) 

Pvac ^universe 

where t un iverse ~ 14 Gyr is the age of the universe. 

If this limit is saturated then we have a model for the ultra-high energy cosmic rays. 
The spectrum in (5.2) implies a hard energy spectrum dN/dE ~ 10 _Q with a ~ 1.2 at the 
relevant energies. The observed flux spectrum dJ/dE is modified by the rapid decrease of 
the attenuation length by almost a factor of 100 between 4 x 10 10 < E < 10 11 GeV (the 
GZK cutoff). The frequently plotted E 3 dJ/dE thus rises for E < 4 x 10 10 GeV, drops for 
4 x 10 10 < E < 10 11 GeV, and continues to rise for E > 10 11 GeV. This behavior is a 
simple consequence of a hard initial spectrum of protons produced uniformly throughout 
space, and is quite consistent with the data above 10 10 GeV. Below 10 10 GeV the observed 
flux rises much faster with decreasing energy, and some other source for cosmic rays must 
dominate. 
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This picture is similar to models involving decay of supermassive long-lived particles 



||13|| |fl4f ||15|| . The primary difference is that the supermassive particles tend to congregate 
in the galactic halo, and this produces a galactic component to the signal in addition to 
a possible extragalactic one. But the galactic component would not produce a dip in the 
spectrum, and it would produce some amount of large scale anisotropy []14| flUR . In our 



picture the absence of a galactic component is natural, although there remains the question 
of just how spatially uniform the decay mechanism would be. 

The production of cosmic gluon jets also ties in with another mechanism for producing 
air showers above the GZK cutoff, this time involving jets produced outside the £ 3 volume. 
This is because gluon jets contain neutrinos, and these energetic cosmic neutrinos have 
some probability for producing Z-bursts within the £ 3 volume, as they travel toward us |16[] . 
In particular neutrinos with energy within 5E/Er = Tz/Mz ~ 3% of the Z-resonance 
energy = 4(eV/m„) x 10 12 GeV may annihilate with an enhanced cross section on 
the nonrelativistic relic antineutrinos (and vice versa) to produce the Z. We estimate 
(with the same uncertainties as before) that the number of neutrinos in this energy band, 
produced per gluon jet, ranges from 300 to 1100 as Er ranges from 10 13 to 10 11 GeV. The 
probability for such a neutrino to produce a Z-burst within the £ 3 volume is in the range 



0.025% to 1%, depending on the relic neutrino clustering JT^]. On the other hand there 



is an enhancement factor of about 100 for the neutrino flux relative to the direct proton 
flux since the former originates from the whole Hubble volume. Each Z-burst produces a 
couple of protons with typical energies E p « .Er/30, and thus this mechanism produces 
protons peaked in a fairly narrow energy range. It may even be of interest if this peak was 
somewhat below the GZK bound, where the protons-from-jets mechanism was deficient. 



6. Conclusions. 

A warped, compact background space-time with a compact scalar field, introduced to 
address the cosmological constant problem, does not present any obstructions to including 
chiral fermions. Since the compactness of the extra dimension resulted from a periodic 
solution to the field equations, rather than adding brane boundaries, the most appropriate 
approach is to give the universe an S 1 /^ orbifold structure in this dimension. 

A Kaluza-Klein expansion of massless fermion and gauge fields in a warped and os- 
cillating background reveals that the zero mass modes of the fermion do not couple at 
leading order with the oscillating terms in the metric. Thus the gauge fields would most 
readily facilitate the decay of the oscillating gravitational field. 

The decay of the oscillations in the metric should proceed rapidly in the early universe, 
at least until the amplitude is sufficiently small that the oscillatory effects are comparable 
to any de Sitter expansion or to effects from any matter and radiation present, in which 
case the simple picture developed in section 2 breaks down and the subsequent evolution 
is more complex. Yet it would be useful to determine when such oscillations could coexist 
with a familiar cosmology in the 3 + 1 dimensional effective theory since we have seen that 
a small, slowly relaxing oscillation today could provide a source for the observe ultra high 
energy cosmic rays. 
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